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Abstract 

In 1999 Berry and Keating showed that a regularization of the ID classical Hamiltonian H = 
xp gives semiclassically the smooth counting function of the Riemann zeros. In this paper we 
first generalize this result by considering a phase space delimited by two boundary functions in 
position and momenta, which induce a fluctuation term in the counting of energy levels. We 
next quantize the xp Hamiltonian, adding an interaction term that depends on two wave functions 
associated to the classical boundaries in phase space. The general model is solved exactly, obtaining 
a continuum spectrum with discrete bound states embbeded in it. We find the boundary wave 
functions, associated to the Berry-Keating regularization, for which the average Riemann zeros 
become resonances. A spectral realization of the Riemann zeros is achieved exploiting the symmetry 
of the model under the exchange of position and momenta which is related to the duality symmetry 
of the zeta function. The boundary wave functions, giving rise to the Riemann zeros, are found 
using the Riemann-Siegel formula of the zeta function. Other Dirichlet L-functions are shown to 
find a natural realization in the model. 

PACS numbers: 02.10.De, 05.45.Mt, ll.lO.Hi 
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I. INTRODUCTION 



At the beginning of the XX century Polya and Hilbert made the bold conjecture that 
the imaginary part of the Riemann zeros could be the oscillation frequencies of a physical 
system. If true this suggestion would imply a proof of the celebrated Riemann hypothesis 
(RH). The importance of this conjecture lies in its connection with the prime numbers. If 
the RH is true then the statistical distribution of the primes will be constrained in the most 



favorable way 



, |2|. Otherwise, in the words of Bombieri, the failure of t 



(see also {4, [s 



le RH would create 
for reviews on the 



havoc in the distribution of the prime numbers 
RH). 

After the advent of Quantum Mechanics, the Polya-Hilbert conjecture was formulated as 
the existence of a self-adjoint operator whose spectrum contains the imaginary part of the 
Riemann zeros. This conjecture was for a long time regarded as a wild speculation until the 
works of Selberg in the 50's and those of Montgomery in the 70's. Selberg found a remarkable 
duality between the length of geodesies on a Riemann surface and the eigenvalues of the 
Laplacian operator defined on it [of. This duality is encapsulated in the so called Selberg 
trace formula, which has a strong similarity with the Riemann explicit formula relating the 
zeros and the prime numbers. The Riemann zeros would correspond to the eigenvalues, and 
the primes to the geodesies. This classical versus quantum version of the primes and the 
zeros is also at the heart of the so called Quantum Chaos approach to the RH. 

Quite independently of Selbergs work, Montgomery showed that the Riemann zeros are 
distributed randomly and obeying locally the statistical law of the Random Matrix Theory 



(RMT) 



10|] . The RMT was originally proposed to explain the chaotic behaviour of the spec- 



tra of nuclei but it has applications in another branches of Physics, specially in Condensed 
Matter [l^ . There are several universality classes of random matrices, and it turns out that 
the one related to the Riemann zeros is the gaussian unitary ensemble (CUE) associated to 
random hermitean matrices. Montgomery analytical results found an impressive numerical 
confirmation in the works of Odlyzko in the 80's, so that the CUE law, as applied to the 
Riemann zeros is nowadays called the Montgomery-Odlyzko law 12|. An important hint 
suggested by this law is that the Polya-Hilbert Hamiltonian H must break the time reversal 
symmetry. The reason being that the CUE statistics describes random Hamiltonians where 
this symmetry is broken. A simple example is provided by materials with impurities subject 



2 



to an external magnetic field, as in the Quantum Hall effect. 

A further step in the Polya-Hilbert-Montgomery-Odlyzko pathway was taken by Berry 
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14| . who noticed a similarity between the formula yielding the fluctuations of the number 



of zeros, around its average position En ~ 27m/ logra, and a formula giving the fluctuations 
of the energy levels of a Hamiltonian obtained by the quantization of a classical chaotic 
system 15|. The comparison between these two formulas suggests that the prime numbers 



p correspond to the isolated periodic orbits whose period is log p. In the Quantum Chaos 
scenario the prime numbers appear as classical objects, while the Riemann zeros are quantal. 
This classical/quantum interpretation of the primes/zeros is certainly reminiscent of the one 
underlying the Selberg trace formula mentioned earlier. A success of the Quantum Chaos 
approach is that it explains the deviations from the CUE law of the zeros found numerically 
by Odlyzko. The similarity between the fluctuation formulas described above, while rather 
appealing, has a serious drawback observed by Connes which has to do with an overall sign 
difference between them [l^. It is as if the periodic orbits were missing in the underlying 
classical chaotic dynamics, a fact that is difficult to understand physically. This and other 
observations lead Connes to propose an abstract approach to the RH based on discrete 



mathematical objects known as adeles 



16|. The final outcome of Connes work is a trace 



formula whose proof, not yet found, amounts to that of a generalized version of the RH. 
In Connes approach there is an operator, which plays the role of the Hamiltonian, whose 
spectrum is a continuum with missing spectral lines corresponding to the Riemann zeros. 
We are thus confronted with two possible physical realizations of the Riemann zeros, either 
as point like spectra or as missing spectra in a continuum. Later on we shall see that 
both pictures can be reconciled in a QM model having a discrete spectra embedded in a 
continuum. 

The next step within the Polya-Hilbert framework came in 1999 when Berry and Keating 



17 



on one hand and Connes [l6l] on the other, proposed that the classical Hamiltonian 
H = xp, where x and p are the position and momenta of a ID particle, is closely related 
to the Riemann zeros. This striking suggestion was based on a semiclassical analysis of 
H = xp, which led these authors to reach quite opposite conclusions regarding the possible 
spectral interpretation of the Riemann zeros. The origin of the disagreement is due to the 
choice of different regularizations of H = xp. Berry and Keating choosed a Planck cell 
regularization in which case the smooth part of the Riemann zeros appears semiclassically 
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as discrete energy levels. Connes, on the other hand choosed an upper cutoff for the position 
and momenta which gives semiclassically a continuum spectrum where the smooth zeros are 
missing. All these semiclassical results are heuristic and lack so far of a consistent quantum 
version. It is the aim of this paper to provide such a quantum version in the hope that it 
will sed new light concerning the spectral realization of the Riemann zeros. 

The organization of the paper is as follows. In section II we review the semiclassical 
approaches to H = xp due to Berry, Keating and Connes which give an heuristic derivation 
of the asymptotic behaviour of the smooth part of the Riemann zeros. Then, we generalize 
the semiclassical Berry-Keating Planck cell regularization of xp by means of two classical 
functions which define a wiggly boundary for the allowed semiclassical region in phase 
space. This generalization allow us to explain semiclassically the fluctuation term in the 
spectrum. In section III we define the quantum Hamiltonian associated to the semiclassical 
approach introduced above. The Hamiltonian is given by the quantization of H = xp 
plus an interaction term that depends on two generic boundary wave functions associated 
to the classical boundary functions of the semiclassical approach. In section IV we solve 
the Schroedinger equation finding the exact eigenf unctions and eigenenergies in terms of a 
function J-'{E) which plays the role of a Jost function for this model, and whose analyticity 
properties are studied in section V. In section VI we find the boundary wave functions that 
give rise to the quantum version of the semiclassical Berry-Keating model for the smooth 
zeros of the Riemann zeta function, which are common to all the even Dirichlet L-functions. 
We also find the boundary wave functions associated to the smooth approximation of the 
zeros of the odd Dirichlet L-functions. In section VII we quantize the relation between 
the fluctuation part of the spectrum and the semiclassical phase boundaries, obtaining the 
equations satisfied by the boundary wave functions, and we solve them explicitely. Finally, 
using the duality properties of these wave functions and the Riemann- Siegel formula of the 
zeta function we find a model whose Jost function is proportional to the zeta function. From 
this fact, and making some additional asumptions, we show that the Riemann zeros on the 
critical line are bound states of the model. However we cannot exclude the existence of zeros 
outside the critical line, which would imply a proof of the RH. We describe in an appendix 
the computation of the wave functions associated to the smooth and exact Riemann zeros. 



The present work is closely related to those in references 



19 



2G 



2l|, where we studied 



an interacting version of the xp Hamiltonian based on the relation of this model with the so 
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called Russian doll model of superconductivity 



22 



23 



to the RH inspired by the latter works see reference 



24] . For a field theoretical approach 
25| . We would like also to mention 



some important differences between the present paper and those of references 



19 



19 



20 



2l|. 



20,1211 to belong to the finite interval 



First of all, the position variable x was choosen in 
(1, N) with oo, while in this paper we choose the half line (0, oo) which gives a more 

symmetric treatment between the position and momentum variables. Secondly, in the earlier 
references the interaction term was added to the inverse Hamiltonian l/{xp), while in this 
paper we add the interaction directly to the Hamiltonian xp, which is more natural from a 
physical viewpoint. We have also tried to make an extensive use of the duality symmetry of 
the Riemann zeta function refiected in the functional relation it satisfies. 



II. SEMICLASSICAL APPROACH 



The classical Berry-Keating-Connes (BKC) Hamiltonian [16|, ll7|, ll8fl 

H^^ = xp, (2.1) 

has classical trayectories given by the hyperbolas (see fig. la) 

x{t)=xoe\ p{t)=poe-'. (2.2) 

The dynamics is unbounded, so one should not expect a discrete spectrum even at the 
semiclassical level. To overcome this difficulty. Berry and Keating proposed in 1999 to 
restrict the phase space of the xp model to those points (x,p) where |x| > Ix and \p\ > Ip, 
with Ixlp = 2Tih. These constraints lead to a finite number of semiclassical states, Af{E), 
with energy between and E given by 

where A is the area of the allowed phase space region below the curve E = xp. The result, 
in units h = 1, is 

A/-B.W = ^(log^-l)+l (2.4) 

which agrees with the asymptotic limit of the smooth part of the formula giving the number 
of Riemann zeros whose imaginary part lies in the interval (0, E), 

imE))~^L,^^-l)+l^OiE-^ (2.5) 




FIG. 1: la) a classical trayectory (j2.2p . The region in shadow is the allowed phase space of the 
semiclassical regularizations of Berry and Keating, lb) generalization of the phase space region 
given by equations (|2.12|) 

The exact formula for the number of zeros, A/ij(-E), due to Riemann, also contains a fluctu- 
ation term which depends on the zeta function |l| (see figl2]), 

^^n{E) = {^^iE))+^^^iE) (2.6) 
mE)) = ^ + 1 

7r 

AAfl(E) = ^ImlogcQ + ^i?) 
where 0{E) is the phase of the Riemann zeta function ({1/2 — iE), 

9{E) = Im log r Q + - I log TT (2.7) 
whose asymptotic expansion 

E, f E\ E TT 



yields (12.51) . The function C,{s), for Re s > 1, can be related to the prime numbers p thanks 
to the Euler product formula 

P>1 ^ 
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FIG. 2: Number of Riemann zeros in the interval (0,-E'): black: exact formula (j2.6p . red: smooth 
function {M{E)), blue: {M{E)) + 1/2. 

This expression diverges if Re s = 1/2, however one can heuristically use it to write the 
fluctuation term in (12. 6p as 

1 °° 1 

AffiiE) = sin(mElogp) (2.10) 

p m=l 

which gives a reasonable result after truncating the sum over the primes. As observed by 
Berry, eg. (12. 101) resembles formally the fluctuation part of the spectrum of a classical ID 
chaotic Hamiltonian with isolated periodic orbits 

1 °° 1 
WE) = - E E o • u X m ^i^(^ci(i?)) (2.11) 



TT 

7p 



m 2 sinh(mAp/2) 



where 7p denotes the primitive periodic orbits, the label m describes the windings of those 
orbits, ±Ap are the instability exponents and Sc\{E) is the classical action, which is equal 
to mET^^, with T^^ the period of 7p. Comparing (12.101) and (12. lip . Berry conjectured the 
existence of a classical chaotic Hamiltonian whose primitive periodic orbits would be labelled 
by the prime numbers p = 2,3,..., with periods Tp = logp and instability exponents 



A„ = ±logp 



13|, 



14| . Moreover, since each orbit is counted once, the Hamiltonian must 
break time reversal (otherwise there would be a factor 2/7r in front of eq. (I2.10p instead 
of 1/vr). The quantization of this classical chaotic Hamiltonian would likely contain the 



7 



Riemann zeros in its spectrum. This idea is the key of the Quantum Chaos approach to the 
Riemann hypothesis. 

Besides the fact that the earher Hamiltonian has not yet been found there is the Connes 
criticism that the similarity between eqs. fl2.10p and (12.111) fails in two issues. The first is 
the overall minus sign in fl2.10p as compared to (12.111) . and the second is that the term 
2 sinh(mAp/2) only becomes p'"/^ when m oo. Connes relates the minus sign problem to 
an alternative interpretation of the Riemann zeros as missing spectral lines as opposed to 
the conventional one (we shall come back later to these conflicting interpretations). These 
two problems were the main Connes's motivations to develop the adelic approach to the RH. 

As we saw above, the Quantum Chaos approach suggests that the fluctuation part of 
the spectrum of the yet unknown Riemann Hamiltonian has a classical origin related to the 
prime numbers. Taking into account the Berry-Keating heuristic derivation of the smooth 
part of the spectrum, it is tempting to extend the semiclassical approach in order to explain 
the fluctuation term in the Riemann formula for the zeros. The simplest idea is to generalize 
the allowed phase space of the xp Hamiltonian replacing the boundaries = Ix and \p\ = Ip 
by two curves Xci{p) and Pci{x), such that (see fig lb) 

x>Xd{p), \p\>Pciix) (2.12) 

where Xci{p) and Pc\{x), are positive functions satisfying 

Xciip) = Xdi-p) > 0, Vp G ]R (2.13) 
Pci{x) > Vx G R+ 

These conditions split the allowed phase space into two disconnected regions in the first 
and forth quadrants of the xp plane. Notice that x is always positive while p can be either 
positive or negative. The BK boundaries obviously correspond to the choice 

BK : Xci(p) = /x, Pdix) = Ip (2.14) 

For the extended BC's the minimal distance Ix and minimal momentum Ip can be defined as 
the intersection point of the curves, Xci{p) and Pci{x), which we shall assume to be unique, 
and satisfying 

Xcl{lp) = Ix, Pc\{lx) = Ip (2.15) 
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The classical xp Hamiltonian together with the BK conditions have the exchange symmetry 

f - f (2.16) 

Ix hp 

whose generalization to the extended model is 

Xc\{lpx/lx) Pel (a;) 



Ix Ip 



(2.17) 



The counting of semiclassical states is based again on eq. (12.31) . The area below the curve 
E = xp and bounded by the conditions (12.121) is given by (see fig. lb) 

l>xj rlpx/lx rXM fE/x 

A = dx dp+ dx dp (2.18) 

''Ix JPclix) Jxi JPcl{x) 

cpi rlxp/lp rPM pE/p 



+ dp dx + dp dx 

Jlp Jxcl{p) Jpi Jxclip) 

The quantities xm, Pm (resp. Xj,pi) are the position and momenta of the points where the 
curve E = xp intersects the boundaries Pci{x), Xci{p) (resp. the line x/l^ = p/lp), and satisfy, 

X J Pj 

E = xm Pci{xm) = Xci{pm) Pm = xipi, — = — (2.19) 

Ix Ip 

The integration of (I2.18P yields 

A = E\og (^-^^ + E- Up (2.20) 
- E\og(P-^)-E\oJ^-^ 

\ ''p J \ '•X 

\i rPM 
dxpc\{x) - / dp Xc\{p) 



Partial integrating the last two terms in (I2.20p and dividing hy h = Ixlp = 2TT{h = 1), the 
semiclassical value of N'{E) reads 

U{E) = Ubk{E) (2.21) 

f Pci{xm)\ E fXc\{pM] 
log - 7^ log 



2vr \ ^p J 27r \ 

r^' dx dpJx) r^" dp dxJp) 

+ — X — h / — P- 

Jlx 



2tt dx Ji^ 2n dp 

The BK conditions (12.141) of course reproduce eq. (12. 4p . More general boundary functions 
induce a fiuctuation term in the counting formula of a form which recalls eq. (l2.6p . Let us 
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denote this term as 



r^' dx dpa{x) dp dxM 

Ju 27r ^ dx J, 2tx ^ dp 



so that 



J\r{E)=J\rBK{E)+n^{E) (2.23) 
Taking the derivative of f l2.22p with respect to E, and using eqs.f l2.19]) one gets 

dn,{E) _ 1 / Pci(xm)\ _ ^ log / ^I^m) ^ (2.24) 



rf^ 27r \ L / 27r \ / 



which imphes that the boundary functions are related to the fluctuation part of the density 
of states. A further simplification is achieved imposing the xp symmetry (12.171) 

Pc\{.xm) _ Xc\{pm) Pm _ xm ^2 25) 



/ 111 



dn,{E) 1 / p_^\ ^ _ 1 / 0^ ^2.26) 



which leads to 



dE TT \ J \ Ix 

Hence, xp-symmetric boundary functions Pci{xm) and Xci{pm) are completely fixed by the 
density of the fluctuations. To find Pci{x), one combines (12.261) and (I2.19P 

Pc,(xm) = e-'^-n^^) = ^, nm = (2.27) 

Xm dE 

which gives xm as a function of E 

XM = y e""fl(^) (2.28) 

Lp 

If nf[{E) = , the latter equations reproduce the BK boundary conditions (12.151) . Eq. (l2.28p 
gives Xm as a function of E and it is monotonically increasing provided 

dxMiE) d'^naiE) 

^>0^1+.i5-Jji>0 (2.29) 

Under this condition we can expressed as a function of xm and replaced it in (12.270 . 
obtaining the boundary function p^^^ = E{xm)/xm- In this case the inverse problem of 
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finding a Hamiltonian given the spectrum has a unique solution at the semiclassical leveL If 
the fluctuations are strong enough at some energies, then condition (12.291) could be violated 
implying that E = E{x) as well as Pc\{x) will be multivalued functions. This gives rise to a 
manifold of boundary functions, each one having discontinuities at some values of x. 



III. FROM CLASSICAL TO QUANTUM 



In this section we shall give a quantum version of the semiclassical results obtained above. 
The starting point is the quantization of the classical hamiltonian Hq = xp. Let us consider 
the usual normal ordered expression 



Ho = ^{xp + px) 



d 1 



(3.1) 



where p = —id/dx. In references 20|, |26|] it was shown that Hq becomes a self-adjoint 
operator in two cases where the domain of the x variable are choosen as: 1) < x < oo or 
2) a < X <h with a and h finite. For the purposes of this paper we shall confine to the case 
1. Case 2 was discussed at length in [20|]. Since a; > one can write (13. ip as 



Ho = x"^px^'\ 



X > 



(3.2) 



The exact eigenfunctions of (13. 2p are given by 

0s(x) = 



^ G ]R 



(3.3) 



where the eigenenergies E belong to the real line. The normalization of (13. 3p is the appropiate 
one for a continuum spectra. 



/■oo 

{(I>e\<I>e') = / dx 0^(x)0£;.(x) = 5{E - E'). 
Jo 

The quantum Hamiltonian associated to the semiclassical approach is 



(3.4) 



(3.5) 



where ipa and ipb are two wave functions associated to the boundary functions Pc\{x) and 
Xci{p), respectively, i.e. 

'4'a ^ Pclix), ^ Xc\{p) (3.6) 
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(3a) 



(3b) 



FIG. 3: Graphical representation of the classical transport operation in the phase space of the xp 
model given in eqs. (|3.9p and (jS.lSp . 

We shall choose real functions il>a{x) and ipbi^) so that H is an hermitean and antisymmetric 
operator, which implies that the eigenvalues appear in pairs E, ~E. The interaction term in 
(13. 5p can be justified by the following heuristic argument. Let us consider a particle which 
at t = belong to the classical allowed region, i.e. xq > Xci(po) and po > Pci{xo). According 
to the classical evolution (12. 2p . the position x{t) increases while the momenta p{t) decreases, 
i.e. 

Classical evolution ■.(xq,pq) — ^(e*a;o, e~*po) (3.7) 
until a time tM where the particle hits the Pci-boundary. 

(e**^Xo, e-*^Po) = {xm,Pci{xm)) (3.8) 

The semiclassical approach suggests to transport this particle from the Pci-boundary to a 
point in the Xd-boundary, (see fig. [3]) 

Classical transport ■.{xm,Pci{xm)) {xci{pm),Pm), (3.9) 

while preserving the energy, 

E = xopo = xm Pc\{xm) = Xc\{pm) Pm (3.10) 

Equation flS.lOp coincides with (12.190 if we choose (xo,Po) = {xj,pi). The transported 
particle at the x^- boundary continues its classical evolution returning to the initial point 
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{xo,po) after a time 

te = — log = — log (3 11) 

This is also the period of the classical trajectory which has become a closed orbit thanks to 
the transport operation f l3.9p . The semiclassical calculation of the previous section measures 
classical action associated to this periodic orbit. At the quantum level the free evolution of 
a state ip is given by the unitary transformation 

Quantum evolution #(0)) — >\i;{t)) = e"^*-^" |V^(0)) (3.12) 

The operator that performs the transport (13.91) is given by one of the interacting terms in 
the Hamiltonian (13. 5p . 

Quantum transport :|'?/') —i\ipb){4'a\i') (3.13) 

which consists in the proyection of the state into the quantum state ipa, yielding the state 
ipb as a result. The hermiticity of the Hamiltonian H implies the existence of the inverse of 
the process (I3.13p . i.e. 

Inverse quantum transport -.lip) —>■ —i\ipa){'4'b\'4') (3-14) 

whose classical analogue is (see fig. [3]d), 

Classical Inverse transport : {xci{pm),Pm), {xm,Pci{xm)) (3.15) 

What is the physical meaning of this process? Let us take for a while a particle in the 
classical forbbiden region where xq < Xci{po) but po > Pci{xo)- This particle will evolve 
freely according to eqs. (l3.7p . until a time tjv/ where it hits the Xd-boundary, i.e. 

(e**^Xo, e-*>o) = Mpm),Pm)) (3.16) 

Then one can apply the inverse transport (I3.15P which carries the particle to the Pci-boundary 
where it continues its free and unbounded evolution : x ^ oo and p 0. The phase space 
area traced by this trayectory is infinite which implies that the number of these kind of 
semiclassical states is infinite forming therefore a continuum. 

In summary, the transport operations between the two boundaries leads classically to 
closed periodic trayectories in the allowed phase space and to open trayectories in the forb- 
biden region. Semiclassically the closed periodic trayectories give rise to bound states while 
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the open ones form a continuum. This is scenario that comes out from the solution of the 
quantum model, as we show in the next section. 

The existence of a semiclassical continuum in the xp model was proposed by Connes in 



reference 16j. Instead of the boundary conditions set by and Ip, Connes restricts the 
phase space of the model to be |x| < A, \p\ < A, where A is a cutoff which is sent to infinite 
at the end of the calculation. The number of semiclassical states is given now by 

XciE) = ^\ogA-^ (log ^-l) (3.17) 

where the first term leads, in the limit A — * oo, to a continuum while the second term coin- 
cides with minus the average position of the Riemann zeros (12.41) . A possible interpretation 
of these result is that the Riemann zeros, are missing spectral lines in a continuum, which 
is in apparent contradiction with the Berry-Keating interpretation of the zeros as bound 
states. As we shall show below both interpretations can be reconciled at the quantum level 
where the Riemann zeros appear as discrete spectra embbeded in a continuum of states. 

IV. EXACT SOLUTION OF THE SCHROEDINGER EQUATION 

In this section we shall find explicitely the eigenstates and the eigenergies of the Hamil- 
tonian (13.51) for generic states ipa and ipb. The method used is similar to the one employed 



in reference 



20| , where instead of the Hamiltonian xp we added an interaction to 1/ xp. The 



Schroedinger equation for an eigenstate iPe{x) with energy E is given by 

- i (x-^ + M i^Eix) + i {Mx){Hi^E) - M^){^a\^E) = E^e{x)) (4.1) 



dx 2^ 

Let us introduce the variable q 

g = logx, gGR (4.2) 

and the overlap integrals 

/■oo 

A = {^Pa\^E)= / dx iJaix)ijEix) (4.3) 







oo 



B = {ipbl-ipE) = dx ilJb{x)ipE{x) 
Jo 

which depend on E. Using these definitions eq. (l4.1l) becomes 
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The general solution of this equation is given by 



Mq) = e-(i/2-^)^ 



Co+ / dq' e 



/ n/2~iE)q' 



where Co is an integration constant. It is convenient to define the functions 



a{q) = e^'^^ipaiq), V'a(a;) 

b{q) = e'/'Mq), M^) 



b{x) 



X 



so that 



Co+ r dq' e-'^^^'iB a{q')-Ab{q')) 

J — oo 



An alternative way to express (14. 7p is 



^^(g) = e-(l/2-.i?k 

where Coo is related to Co by 



Cno 



oo 



, dq' e-'^^ {B a{q') - A b{q')) 
'1 



Coo = Co + 5 a{-E) - A b{-E) 



where 

/oo 
dqe'^^^fiq), f = a,b 
-00 

We shall assume that a{q) and 6(g) satisfy 



lim 



r_^dq'e-^^^'f{q') = 0, f = a,b 



lim,^oo /,°° dq' e-^^^'f{q') =0, f = a,b 

which implies that the asymptotic behaviour of 'ipEix) is dominated by Co, Co 

C C 
lim ipEix) = I , lim = — 7^^, 

Plugging fl4.7l) into (14 ■3p yields the relation between the constants A, B, Co, 



Sbfi 1 — Sb.a I \ B 



Co 



a{E) 

m 
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where the functions Sf^g{E) with f^g = a^b are defined by [27| 

/oo pq 
dqe'^^fiq) / dq' e-''"^' giq') 
-oo J —oo 

Similarly, introducing (14.81) into (14 .Sp yields 



(4.14) 



1 Sa,b San 



Sl,,b 1 + Sb,a I \B \ b{E) 



(4.15) 



where 



dq e^^^ /(g) 



This function is related to Sfn in two ways 



(4.16) 

(4.17) 
(4.18) 

To derive these equations one makes a change of order in the integration. Combining (I4.17p 
and (I4.18P one obtains the shuffle relation 



Sf,,{E) = -S^E) + f{E) g{-E) 
SUE) = S,j{-E) 



SUE) + S,j{-E) = f{E) g{-E) 



(4.19) 



The terminology is borrowed from the theory of multiple zeta functions where there is 
a similar relation between the two variable Euler-Zagier zeta function (^(si,S2), and the 
Riemann zeta function C{s) 



28 



29|. 



The solutions of the eqs. (l4.13p and fl4.15p depend on the determinant of the associated 
2x2 matrices given by 



T{E) — 1 + Sa,b — Sb,a + Sa,aSb,b — Sa,bSb,a 
T{E) = I — Sa,b + Sb^a + Sa,aSb,b — Sa^bSb,a 



(4.20) 



which are related by (14.181) 



J^{E) = J^i-E) 



Moreover, since a{x) and b{x) are real functions one has 



which in turn implies 



SUE) = SU-El 

T*{E) = J^{-E*) 
16 



(4.21) 



(4.22) 



(4.23) 



After these observations we can return to the solution of fl4.13p and fl4.15p . We shall distin- 
guish two cases: 1) J-'{E) ^ and 2) T{E^ = 0, where E is real since it is an eigenvalue of 
the Hamiltonian fl3.5p . 
Case 1: J^{E) ^ 

Eq. fl4.23p implies that J-'{—E) ^ and therefore A and B can be expressed in two different 



ways, 



A 



B 



HE) 

Coo 

H-E) 

Co 
HE) 

Cna 



1 - S,^a) a{E) + Sa,a h{E) 
(1 + 5,,,) a{E)-Sa,ah{E) 
Sb,b a{E) + (1 + Sa,b) b{E) 
Sb,b a{E) + {l-Sa,b)b{E) 



H-E) 

Now using eq. (l4.17p . these eqs. reduce to 

Co ^(^) 



(4.24) 



(4.25) 



(4.26) 



Coo H-E) 

which by eq. fl4.23p is a pure phase for E real. Hence, up to an overall factor, the integration 
constants for this solution can be choosen as 

Co = J'iE) 

Coo = H-E) (4.27) 

A = {l-Sb,a)a{E) + Sa,ab{E) 

B = -Sb,ba{E) + {l + S,^b)b{E). 

Since the constants Cq, Coo do not vanish, the wave function is non normalizable near the 
origin and infinity (recall eq. fl4.12p ) and therefore they correspond to scattering states. Of 
course they will be normalizable in the distributional sense. 
Case 2: J^{E) = 0. 

The integration constants can be choosen as 



Co 

Coo 

A 
B 






Sa,a 

(1 + Sa,b) 



(4.28) 
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Bound states 



Scattering states 




FIG. 4: Pictorial representation of the spectrum of the model. The bound states are the points 
where J^{E) = 0, which are embedded in a continuum of scattering states. 

which solves eqs. (14.131) and (I4.15p . Since Co = Coo = 0, the leading term of the behaviour 
of iPe{x) vanish near the origin and infinity and under appropiate conditions on ipa,b, the 
state ipE will be normalizable corresponding to a bound state. In the appendix we compute 
the norm of these states. 

Hence the generic spectrum of the Hamiltonian (13.51) consist of a continuum covering 
the whole real line with, eventually, some isolated bound states embedded in it, whenever 



J^{E) = 0. This structure also arises in the Hamiltonian studied in reference [20||. The 
function J-'{E) plays the role of the Jost function since its zeros gives the position of the 
bound states and its phase gives the scattering phase shift according to eg. (14.261) . 

Before we continue with the general formalism it is worth to study a simple case which 
illustrates the results obtained so far. 

An example: a quantum trap 

Let us start with the classical version of a trap where a particle is restricted to the region 
Xb < X < Xa- The semiclassical number of states is given by the area formula (12.31) . 

A dx E E ^ Xa , , , 

n = — = = _log— (4.29 

2n J^^ 2% X 2n Xb 

which yields the eigenenergies 

^ 2n n , ^ 

En = , , , new. (4.30 

\0g{Xa/Xb) 

The quantum version of this model is realized by two boundary states 4'a,b{x) proportional 
to delta functions, i.e. 

ijja{x) = ao xl^'^6{x - Xa), ^b{x) = K xl^'^6{x - Xb). (4.31) 
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I 




FIG. 5: Semiclassical picture of the model represented by the potential (j4.32p . 

The associated potentials a{q),b{q) are 

a{q) = ao5(g - Qa), b{q) = bo6{q - Qb), (4.32) 

Qa = \0gXa, Qb = log Xfc 

The various quantities defined above are readily computed obtaining 



(4.33) 



a 


= ao e'^^\ 


b= 


Sa,a 


- ^b,b 


_ 

2 


Sa,b 


= ao bo e*^^-" 


Si 


where qa,b = Qa - Qb = log{xa/xb). 


Plugging these 


eqs. 


HE) 




+ aol 



(4.34) 



For generic values of ao,&0) the Jost function fl4.34p never vanishes obtaining a spectrum 
which is continuous. However, J^{E) vanishes provided the following condition holds 

e = ^ = ±1 ^ jr(^) = 2(1 + e e'^ ^'"■*) (4.35) 

in which cases the spectrum contains bound states embbeded in the continuum with energies 

Ife= 1^^ = '"^^ + ^/'^ nGN (4.36) 

Cla,b 

^ 2™ 
If e = -1 ^ E„ = n e IN 

Qa,b 
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that agree with the semi classical energies fl4.30p for n » 1. The unnormalized wave function 
of the bound states, i.e. J-'{E) = 0, can be computed from eq. (14. 7p 



1 \ 1, Xb < X < Xa 

M^) = -:tj^-{' (4.37) 

' U, X < Xb or X > Xa 



X 



which shows that they are confined to the region {xb, Xa)- The wave functions when J^{E) ^ 
can be similarly found obtaining 



J^{E), < X < Xb 

Xb<X<Xa (4.38) 

J-'{ — E), Xa < X < oo 



Hence if (I4.35P holds, these wave functions vanishes in the region [xb-, Xa) which contains 
the trapped particles (14.370 . In this example the mechanism responsible for the existence 
of bound states is the transport of the particles from the position Xa to the position Xb- At 
the quantum level the confinement requires the fine tuning of the couplings (see eq. (14.351) ). 
which introduces periodic or antiperiodic boundary conditions depending on the sign of e. 
When |e| 7^ 1 the particle can scape the trap and the bound states become resonances. 



V. ANALYTICITY PROPERTIES OF F{E) 



As in ordinary Quantum Mechanics, the Jost function J^{E) satisfy certain analyticity 
properties reflecting the causal structure of the dynamics. In our case these properties 
follows from those of the function Sj^g (eq. (I4.14p ) and the definition (I4.20p . 

Indeed, let us express Sf^g{E) in terms of the Fourier transforms of the functions f,g. 
First we replace g{q) by its inverse Fourier transform 



back into eq. (l4.14p . obtaining 



dE', 



g{-E') / dq e^^V(g) / dq' e^^^'"^)^' 



The last integral is given by the distribution 



dq e 



I ^i(E'-E)q' ^ ^iq{E'-E) 



rc6{E'-E) + -P-^ 



E' -E 



(5.1) 



(5.2) 



(5.3) 
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where P denotes the Cauchy principal part. Plugging (15.31) into (15. 2p and using the Fourier 
transform of / gives, 



f{E) g{~E)+P r 

J — c 

Alternatively, one can write (15. 4p as 



dE' fjE') g{-E') 
m E' -E 



(5.4) 



dE' fjE') g{-E') 



oo 



with e > an infinitesimal. Eq. (15.51) shows that the poles of Sf^g{E) are located in the 
lower half of the complex energy plane. Thus for well behave functions /,^, the function 
Sf^g{E) will be analytic in the complex upper-half plane. These properties also apply to 
J-'{E) which is the product of Sf^g functions with f,g = a, b. Another important property 
of the Jost function T{E) is that its zeros lie either on the real axis or below it, i.e. 

If T{E) = ^ Im E < (5.6) 



The proof of this equation is similar to the one done in reference 20[, being convenient to 
regularize the interval x G (0, oo) as (A^~^, A^) with A^ — > oo. 

In the appendix we use the results obtained in this section to compute the norm of the 
eigenstates. 

VI. THE QUANTUM VERSION OF THE BERRY-KEATING MODEL 

Let us consider the BK constraints x > and \p\ > Ip. It is rather natural to associate 
constraint x > Ix with the wave function 

V'b(x) = 6o 5(x - /.) (6.1) 

1/2 

which is localized at the boundary x = Ix- The factor IJ gives the correct dimensionality 
to ipb{x), with bf) a dimensionaless parameter. On the other hand the constraint \p\ > Ip 
admits two possible quantum versions, 

=2a.(M"%/'=°*^' (6.2) 

ipg^ (x) V^TT/ I sin(/pa;) 



21 



Due to the fact that ipa has to be real, one cannot choose a pure plane wave e*'p^. The 
boundary wave functions (16.11) and (16. 2p are the cosine and sine Fourier transform of each 
other, namely 

V^+(x) 2ao f Ip V^^ f°° ^ , / X I cos{lpXy/l^, 

dy ifjbiy) X < (6.3) 



i'ai^) ^° V2vr/a;y Jo I sin(/pxy//; 

Indeed, extending the domain of ipb{x) according to the parity of ip2 iv = i) o^i^ gets 



bo 



which are the quantum analogue of the classical equations (12.171) . Later on, we shall consider 
more general wave functions il'a,b to account for the fluctuations in the Riemann formula, 
imposing again eq. (16.31) . The relation (16.31) between ip"^ and tpf, must imply a close link 
between their Mellin transforms a±{E) and b{E). To derive it, let us write 

^ r X 2ao / Ip V^' Ha -^me Ha , i ^ \ ^os{lpXy/Q 

a±{E) = / X ^/'+*^ ^^{x) = — [-^\ / dxx ^/'+*^ / dy My)><S 

Jo bo \2ttI^J Jq Jq I sm{lpXy/Q 

(6.5) 

The basic integrals one needs are 



d..-4«-xr"''"' =i(?^r xT. (6-6) 




sin(px) 2 V \p\ 
where 

{-iE T{l/i+iE/2) _ _|_ 

^ r(i/4-i£;/2) ' V-+ , . 

" r(3/4-i£;/2) ' 

The function 0+{E) coincides with the phase of the Riemann zeta function (12. 7p . and more 
generally of the even Dirichlet L-functions, while 9-{E) is the phase factor of the odd 
Dirichlet L-functions. These phases appear in the functional relation of even and odd L 
functions, and they arise in our context from the two possible relations between the boundary 
functions iIj^ and Plugging eq. (l6.6p into (16.50 yields 

where the integral is nothing but b{—E), thus 

■^^^E) = ^(^Y e'^'^^'^^bi-E) (6.9) 



bo \ Ip 
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This important equation reflects the relation (16.31) which in turn is the quantum version 
of the xp symmetry between boundaries. In the BK case, the Mellin transforms of the 
associated wave functions fl6.1l) and fl6.2p are 

ME) = ao (^^ye"''^^^\ m = bo if (6.10) 

which are pure phases, up to overall constants. The Sj^g functions can be readily computed 
using eq. (15.41) . To do so, we first consider the products 

a±{E) a±{—E) = a^, 
b{E)b{-E) = bl (6.11) 
a±{E) b{-E) =aobo e^^^i^^) 
b{E) a±{-E) = aobo e-2^^±(^) 

where we used l^lp = 27r and that 9±{—E) = —6j-{E). The diagonal terms of Sf^g are given 
simply by 



Sa,,a,{E) = ^, S,AE) = j (6.12) 
since the Hilbert transform of a constant is zero, i.e. 

r°° dt 1 

P - =0, EeH (6.13) 

mt- E 

The computation of Sa^^fi and Sb,a± uses the analytic properties of e^*^±^^). Let us focus on 

the case of e^^^^^^^ = e'^^^^^\ This function converges rapidly to zero as \E\ oo in the 

upper half plane, and it has poles at En = i{2n + 1/2) (ra = 0, 1, . . . ) where it behaves like 

(-1)" 2(27r)2" 1 

{2n)\ 2n + 1/2 + iE ^ ' 

We can split e^*^'-'^-' into the sum 



^2ie(E) 



n+{E) + n^{E) (6.15) 



^ ^ (-1)" 2(2vr)^" 1 

^ ' ^ {2n)\ 2n + 1/2 + iE 



n=0 

where Q^{E) is analytic in the upper half plane and goes to zero at +ioo, while Q_^{E) has 
poles in the upper half plane and behaves as 1/E at infinity. The function (E) can also 
be written as 



.{E) = 2 [ dx x-^^^+'^ cos{2nx) (6.16) 
Jo 



r^A -E 1 5 .E 



l + 2iE '4 2 2 4 2 
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where 1F2 is a hypergeometric function of the type (1,2). From the analyticity properties 
of Q± one gets inmediately their Hilbert transform 



00 



PI -^ = ±n^iE), EeU (6.17) 
Tilt — hi 



00 



Hence Saj^,b = Sa,b, as given by eg. (15.41) . becomes 

'° [n^E) + n^E) + n+{E) - n^E)] 



(6.18) 



2 

Similarly one finds 

Sb,a{E) = aobon4-E) (6.19) 

Notice that both functions are analytic in the upper half plane. The Jost function finally 

reads 

2 



J^{E) = l + aoboin+{E)-n4~E)) + 



2 

-{aobofn+{E)n4-E) (6.20) 



In the asymptotic limit \E\ » 1 



n_iE) ~ 1 ^ n^iE) = e2^^(^) + 0(1) (6.21) 



which implies 

HE) = 1 + ao&oe^^^(^) + (^) ' + 0(1) (6.22) 
This Jost function has zeros on the real axis, up to order 1/E, provided 

g = ^ = ±1 ^ ^(^) = 2(1 + e e2^^(^)) + 0(1) (6.23) 
2 E 

The choice e = — 1 reproduces the smooth part of the Riemann formula (12. 6p since, 

e = -1 ^ 1 - e2^^(^) = 1 - e2-<-^(^» = (6.24) 

where E is the average position of the zeros. On the other hand the choice e = 1 leads to 

e = 1 ^ 1 + e2^^(^) = ^ cos e{E) = (6.25) 
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so that the number of zeros in the interval (0, E) is given by 




(6.26) 



which gives a better numerical approximation than the term {Af{E)) that appears in the 
exact Riemann formula fl2.6p (see also fig. 2). In the case of the sine boundary function (16. 2p 
one similarly obtains the smooth part of the zeros of the odd Dirichlet L-functions. 

In summary, we have shown that the semiclassical BK boundary conditions have a quan- 
tum counterpart in terms of the boundary wave functions ipa,b, and that the average Riemann 
zeros become asymptotically bound states of the model or more appropiately resonances. 

VII. THE QUANTUM MODEL OF THE RIEMANN ZEROS 

In section II we showed how to incorporate the fluctuations of the energy levels in the 
heuristic xp model by means of the functions Pc\{x) and Xci{p) which define the boundaries 
of the allowed phase space. These functions are given by eg. (12.261) in terms of the density 
of the fluctuation part of the energy levels. In the quantum model the functions Pci{x) and 
Xci{p) are represented by the wave functions ipa and ipb- Hence it is natural to impose the 
following conditions 



where n'f^{E) = dn^{E)/dE and Hq is the no interacting Hamiltonian (13. ip . The hat over 
X and p stress the fact that they are operators. Eqs. (l7.ip and (17. 2p can be taken as the 
definition of the boundary wave functions. To solve these eqs. let us write them as 




(7.1) 



(7.2) 



{\og\p\ + Xp + Tin'^{Ho))\i^a)= 
{\ogx + K + 71 n'^{HQ))\i)i,) = 



(7.3) 



(7.4) 



Xp = -log/p, = - \ogl, 



x 



(7.5) 



It is convenient to expand the states \ipafi) in the basis ()3.3p 




(7.6) 
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Let us first consider eq. fl7.4p which in the basis (17. 6p becomes 

dE' {<Pe\ \ogx\(t)E') ME') + (A. + vr n'^{E))ME) = (7.7) 
The matrix elements of the operator log x can be readily computed, 

{(t)E\\ogx\(t)E') =-i5\E'-E) (7.8) 
which replaced in (17.71) and upon integration yields 

^dU^ + (A. + TT n'^{E))UE) = (7.9) 
The solution of (17. 9p is simply 

ME) = Mo e'^^^ ^+-"fl(^)) (7.10) 

where ipbfl is an integration constant. The x-space representation of ipb follows from (I7.10p 
and dUD 

/oo i>oo 1 rp 

dE ME) M^) = Mo / -7^ e^^^^ E+.n,iE))^-l/2+^E ^j -^^^ 

Recalling that ipb{x) = b{x)/^/x one gets 

h{x) = Mo r 4^ ^'^'^ E+.n^(E))^^E ^7_^2) 

7-00 v27r 



Observing that h{x) is related to its Fourier transform h{E), as 

dE 



b{x) = I — b{E) x-'^ (7.13) 
one finally obtains 

h{E) = e-'^^^ E+nna(E)) 

where we assumed that nfi{E) is an odd function of E. If tiq^E) = 0, eg. (17. 140 reproduces 
(Km . i.e. 

nfl(E) = ^ b{E) = V2^ ^b,o = bo if (7.15) 
To simplify the notations we shall write (17.141) as 

6(E) = 6o 4^ e-'""«(^) (7.16) 
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Let us now solve the condition (17. 3p for the wave function ipa- We first need to define the 
operator log \p\ acting in the Hilbert space expanded by the functions (j)E {E G IR). In this 
respect it is worth to remember that the operator p = —id/dx is self-adjomt in the real line 
{—oo, oo) and in the finite intervals (a, b), but not in the half-line (0, oo) [30]. However, the 
operator admits infinitely many self-adjoint extensions in the half-line provide the wave 
functions satisfy the boundary condition 

ijj'{0) = K^{0) (7.17) 

where k G R U oo. We shall confine ourselves to the cases where k = and oo, which 
correspond to the von Neumann and Dirichlet EC's respectively, 

K = ^^'(0) =0, (7.18) 
K = oo ^ ^(0) = 

The corresponding eigenstates of the operator with eigenvalues read 

Xp ^ X I ^ °) (7 19) 

Xp V vr 1^ sin(px) {p > 0) 

These basis are complete in the space of functions defined in (x > 0), i.e. 

dp{x>)r x;{x') = 6{x-x'), x,x'>0, rj = ± (7.20) 







The operator log |p| will be defined as |logp^, and therefore admits the same self-adjoint 
extensions as p^. The analogue of eq. (l7.7p reads now 

dE' log 1^ IM ME') + (Ap + TT n'^{E))ME) = (7.21) 

The matrix elements of log \p\ can be computed introducing the resolution of the identity in 
the basis fOQll . 

/•oo 

(0^1 log li^ 10^,) = / dp hgp {<I>e\x;){x;\<I>e') (7.22) 
Jo 

where the overlap of the eigenstates of p^ and Hq are 



oo 



dx cos(px) 



{xt\<pE) = I - X <^ ' (7.23) 

^ I sin(pa;) 
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These integrals were already computed in eq. (l6.6p . and the result is 



(n \-l/2+iE 



Plugging this eq. into fl7.22l) . and performing the integral gives 

{(t)E\ log 1^ \(t)E') = I 5\E' - E)(27r)*(^'-^) ^ne,[E')-e,{E)) ^^^gS) 
which introduced in (17.211) yields a differential equation whose solution is 

The function ipai^x) reads 

/oo poo J rp 

dE = ^,,0 / ^(27r)-''^e-^(^- E+.n,iE)-,2e,iE)) ^-V^+iE 

-oo J —CO V 27r 

(7.27) 

while 

a^jix) = iPafi r 4i(2vr)-''^e-*(^- E+.n^iE)+2e,iE))^^E (7^28) 
J-oo v2tc 

whose Fourier transform is 

a,{E) = V'a,o(27r)i/2+if^e*(^- E+.n^iE)+2e,iE)) (7_29) 
If there are no fluctuations, eq. (]7.29p reduces to 

,2ier,{E)) 

which coincides with eq. (l6.10p . To simplify notations we shall write (I7.29P as 



n^{E) = =^ a,{E) = v^^,,o (f^) e^^''"^^" (7.30) 



Z^(^E) = ao ^^i^n,iE)+2e,iE)) (7 3^) 

The two solutions (17.161) and (I7.3ip satisfy the duality relation (16. 9p and hence the wave 
functions ^/'a^(x) is the cosine or sine Fourier transform of ipb{x) ( see eq. (16.30 ). 

Having found the boundary wave functions for generic fluctuations we turn into the 
computation of the corresponding Jost function. The basic products of the a and b functions 
needed to find the Sj^g functions are similar to eqs. (l6.1ip . 

a±{E) a±{~E) = a^, 
b{E) b{-E) = bl (7.32) 
a±{E)b{-E) = aobo e2^(e±(^)+™fl(^)) 
b{E) a±{-E) = aobo e-2^(e±(^)+-"fl(£^) 
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The diagonal terms of Sj^g are the same as in eq. fl6.12p . i.e. 

Sa^,a^{E) = ^, M^) = | (7.33) 

while the evaluation of the off-diagonal terms depends on the analytic properties of the 
function e'^'^'^^±(^) where 

n^{E)^^-^ + n4E)) (7.34) 

71 

This definition is strongly reminiscent of the Riemann formula (12. 6p . with n±{E) playing 
the role of N'r{E), and n^{E) that of Affi{E). However, we must keep in mind that N'r{E) 
is a step function while we expect nj-{E) to be a continuous interpolating function between 
the zeros. The value of Sa^^b is given by the integral 



^2TTin±{E) _|_ p 



_^ ni t — E 



(7.35) 



We shall make the asumption that e^'^*"*^^) is an analytic function in the upper half plane 
which goes to zero as \E\ — > oo. In this case the Cauchy integral on the RHS of ( 17.35^ is 
equal to e'^'^^'^±i^) and one finds 

5,^,fe(i?) = ao6oe'™±(^) (7.36) 
Similarly Sb^a± vanishes so that the Jost function reduces to 

T{E) = l + a,h, e2™±(^)+(^^y 



(7.37) 



and under the usual choice 

^ = ±l^^(E) = 2(l + ee2™±(^)) (7.38) 



2 

When nfl = the results of the previous subsection showed that e = 1 gives a better 
numerical estimate to the smooth part of the zeros. In the sequel we shall also make that 
choice which implies that the number of zeros of T{E) in the interval (0, E) is 

Uc^m{E) = KUE) + naiE) = n±{E) + ^ (7.39) 

where N'sia{E) was defined in (16.261) for the particular case of the zeta function ({s), which 
corresponds to nj^{E). Equation (I7.39P agrees asymptotically with the semiclassical formula 
( I2.23p . which confirms the ansatz made for the states tpa and ip^. 
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The connection with the Riemann-Siegel formula 



The next problem is to find the function nfi{E), and therefore AfquiE), which gives the 
exact location of the Riemann zeros. Let us consider the case of the zeta function with the 
following choices of parameters 

V = +, e = 1, ao = bo = V2, 1^ = 1, Ip = 27r (7.40) 

which correspond to the potentials (recall fl7.3ip and fl7.16p ) 

a(t) = e^^^^'W+^fl^*)) = eiim+^^it)) (7-41) 

where we skip a common factor ^/2 and denote n(E) = n^{E). These two functions are 
interchanged under the transformation 

a(t) -> e^^^W a(-t) = h{t) (7.42) 
h{t) e2^^(*) b{~t) = a{t) 

so that their sum is left invariant, 

a{t) +b{t) e^''^'\a{-t) + b{-t)) = a{t) + b{t) (7.43) 

The functional relation satisfied by the zeta function implies 

^(1/2 - ^t) ^ e2^^WC(l/2 + ^t) = C(l/2 - zt) (7.44) 

which suggests to relate a + b and ( as 

ai/2- It) =p{t){a{t) + bit)) (7.45) 

where p{t) is a proportionally factor. Using eqs. (17.421) into f l7.45p yields 

C(l/2 -It) = 2 pit) e^^(*) cos(7rn(t)) (7.46) 

This formula can be compared with the parametrization of the zeta function in terms of the 
Riemann-Siegel zeta function Z{t) and its phase 9{t), 

C(l/2 - It) = Z{t) e'^W (7.47) 
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which leads to, 



Z{t) = 2 p{t) cos(7™(t)) 



(7.48) 



This equation is rather interesting since it imphes that the zeros of cos(7™(t)), which give 
the bound states of the QM model, are also zeros of Z{t), of course if p{t) does not have 
poles at those values. Viceversa, the zeros of Z{t) can be zeros either of cos(7m(t)), or of 
p(t), or both. The latter possibility would be absent if the Rieman zeros are simple, as it is 
expected to be the case. 

A first hint on the structure of the functions p{t) and cos(7™(t)) can be obtained using 
the Riemann-Siegel formula for Z{t), 



u{t) 



Z{t) = 2^n-^'^ cos{e{t) -tlogn) + R{t), v{t) 



n=l 




27r 



(7.49) 



where [x] the integer part of x and R{t) is a reminder of order t Combining the last 
two equations one finds 



Z{t) = 2p{t) [cos^(t) cos(7rnfi(t)) - sin^(t) sin(7rnfl(t))] 
^/ N cosli logn) . ^ , , smitlogn) 



(7.50) 



n=l 



n=l 



n 



u(t) 



n=l 



which suggests the following identifications 

/ \ / / N\ cos(tlogn 
p{t} cos{7ma{t)) ~ 2^ — 

p{t) sin(7™fl(t)) 

that can be combined into 



(7.51) 



y-v sm(r logn) 



n=l 



n 



1/2 



/(t)^p(t)e-«W~5^-^ 



(7.52) 



n=l 



The fluctuation function riait) is then given by the phase of f(t), i.e. 



nfl(t) = -Imlog/(t) 

TT 



(7.53) 



In fig. [6] we plot the values of NqMif) that correspond to the approximate formula (17.521) . 
which shows an excelent agreement with the Riemann formula ( 12. 6p . This is expected from 
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FIG. 6: In black: Mr{E), in red: Mqm{E) in the interval (10,40) 

the fact that the main term of the Riemann-Siegel formula already gives accurate results for 
the lowest Riemann zeros. For higher zeros one has to compute more terms of the reminder 
R{t) depending on the desired accuracy. Observe that NqMif) is a smooth function, except 
for some jumps at higher values of t (not shown in fig. [6]) due to the approximation made, 
unlike Afn{t), which is a step function. 

In fig. [7] we plot the values of (17.531) together with those of the fiuctuation part of the 
Riemann formula (12. 6p . i.e. 

Mn{t) = hm\ogc(^^ + i?j (7.54) 

The jumps in A/fi(t) correspond to the Riemann zeros, while those of nfi(t) correspond, either 
to jumps of the function i/(t) appearing in the Riemann Siegel formula (I7.49p . or to those 
points where the curve f{t) cuts the negative real axis in the complex plane. 

We gave in section II a formal expression of eg. (17.541) in terms of prime numbers, eq. 
(I2.10p . which resembles the fiuctuation part (12. lip of a quantum chaotic system. Eq. (l2.10p 
is based on the Euler product formula (12. 9p which is not valid in the case where s = 1/2 + it, 
since Re s > 1 for convergence of the infinite product. The Euler product formula does not 
apply to the truncated sum (I7.52p . however we shall naively try to establish a relationship. 
Let us denote by Pn the n*'^-prime number, e.g. pi = 2,p2 = 3, etc, and by n(x) the number 
of primes less or equal to x. The sum (17.520 involves all integers up to z/(t), which can be 
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FIG. 7: In black: A/'fl(^), in red: n^{E) in the interval (50,80) 

expressed as products of the first /i(t) prime numbers where 

= n(z/(t)), = inf M < (7.55) 

Using these functions we define a truncated Euler product as 

CE(l/2 + ^t)^n ,_ -1/2-.. (7-56) 

It is easy to see that Ce(1/2 + it) is not equal to f(t), for there are terms in (17.561) which 
do not appear in ( I7.52p . although all the terms appearing in the latter sum also appear in 
the former product. The point is that a numerical comparison of these two functions shows 
a qualitative agreement as depicted in fig. [HI Indeed, the minima and maxima of their 
absolute value are located around the same points, and the same happens for the zeros of 
their arguments. The conclusion we draw from these heuristic considerations is that the 
function f{t) contains some sort of information related to the primes numbers although not 
in the form of an Euler product formula as is the case of (e{^/ "2 + it). It would be interesting 
to investigate the consequences of this results from the point of view of Quantum Chaos. 

The Berry-Keating formula of Z{t) 

The main term of the Riemann-Siegel formula (17.49P is not analytic in t due to the 
discontinuity in the main sum. This problem was solved by Berry and Keating who found 

33 




60 70 



iO 90 100 




FIG. 8: Left: in black: \f{E)\, in red: |Ce(1/2 + iE)\ in the interval (50,100). Right: in black: 
Arg f{E), in red: Arg Ce(1/2 + iE). 



an alternative expression for Z[t) 3l|]. The formula is 

oo 

Z{t) = J2iTn{t)+U-t)) 

where 



(7.57) 



n=l 



ut) = T:{-t) 



J e(t) 



(7.58) 



Pnit) = -- j —e 



2Txi 



C- 



and C- is an integration contour in the lower half plane with Im < — 1/2 that avoids a cut 
starting at the brach point z = —t — i/2. The constant K in (17.581) can be choosen at will 
and it is related to the number of terms of the RS formula that has been smoothed for large 
values of t. Using eq. (l7.57p one can write the zeta function as 



C(l/2 it) - e ^ ^ ^ ^i/2+it ' ^i/2-it 

n=l n=l 



(7.59) 



which can be compared with (I7.45P obtaining 



f{t) = p(t)e'™«(*) 



/ J „l/2+it 



(7.60) 



n=l 



SO that (17.591) can be written as 



C(l/2-^t)=e2^^W /W + /H) 



(7.61) 
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Eg. (17.601) gives an exact expression of f{t), which is in fact a smooth version of (I7.52p . Berry 
and Keating also found a series for Z{t) which improves the RS series. The first term of 
that series corresponds to the following value of the Pn(t) functions 

P^^\t) = ^Erfc fi^v^") (7.62) 



2 VQ(^,^) J 
^{n,t) = \ogn-9\t), Q\K,t) = -it9"{t) 

where Erfc is the complementary error function. Using these formulas one can find a better 
numerical evaluation of the functions Mqui^) and n^{t). 

It is perhaps worth to mention that eg. (17.611) . with the approximate value of f{t) given 
by (I7.52p . is a particular case of the so called aproximate functional relation due to Hardy 
and Littlewood Q] 

as) = J2n~^ + vr^-V2 ^^7;0^Y.^'~' + Oix^^ + 0{\t\'/'-y-') (7.63) 

n<x \ I J n<y 

where s = a + it, \t\ = 27Txy, < a < 1. Recalling that in our model t is the energy E, 
then eguation \t\ = 2TTxy becomes the hyperbola \E\ = xp with p = 2TTy = Ipy so that the 
sums in (I7.63P run over the integer values of the positions and momenta in units of and Ip 
respectively. Eg. (l7.63p also suggests that the case where a 1/2 could be related to the non 
hermitean Hamiltonian Hq = {xp + px)/2 — i{cr — 1/2) whose right (resp. left) eigenfunctions 
are given by l/a;°"~*^ ( resp. l/x^~'^~'^^). 

On more general grounds, we would like to mention two important points. First is that 
one still needs to show that the function n(t), defined in eg. (17.340 . is such that e^'^*"^*^ 
is analytic in the upper-half plane and that it goes to zero as \t\ oo, so that the Jost 
function is indeed given by eg. (17.390 . as we have assumed so far. Second, and related to the 
latter point, is that that the function nfi{t) is well defined provided f{t) does not vanish for 
t real, in which case (I7.6ip reads also 

C(l/2 - zt) = f{-t) (l + e'^'^'^jp^) = fi-t) m (7.64) 

which shows that our construction of a QM model of the Riemann zeros relies on the absence 
of zeros of the function f{t) on the critical line. These zeros were investigated by Bombieri 
long ago in an attempt to improve the existing lower bounds for the number of Riemann 



zeros on the critical line 



32| . In this regard our results give further support, but not a proof. 
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to the RH. As suggested in 20|, l21j that proof would follow if the zeta function ({1/2— it) can 



be reahzed as the Jost function of a QM model of the sort discussed so far, due to its special 
analyticity properties. Eq. fl7.64p gives a partial realization of this idea but the function f{t) 
lacks of a physical interpretation so far. The latter approach is analogue to the ones proposed 
in the past by several authors where the zeta function gives the scattering phase shift of 
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some quantum mechanical model, particularly on the line Re s = 1 

Another important question is: where are the prime numbers in our construction? As 
suggested by the Quantum Chaos scenario, the prime numbers may well be classical ob- 
jects hidden in the quantum model, so the next question is: what is the classical limit of 
the Hamiltonian?. The free part is of course given by xp, but the interacting part is an 
antisymmetric matrix with no obvious classical version. The existence of such a classical 
Hamiltonian may help to answer the prime question but it may also lead to a real physical 
realization of the model. Work along this direction is under progress 38 1. 



Acknowledgment s 

I wish to thank for discussions M. Asorey, M. Berry, L.J. Boya, J. Garcfa-Esteve, J. Keat- 
ing, J.I. Latorre, A. LeClair, J. Links, M.A. Martm-Delgado, G. Mussardo, J. Rodriguez- 
Laguna and P.K. Townsend. This work was supported by the CICYT of Spain under the 
contracts FIS2004-04885. I also acknowledge ESF Science Programme INSTANS 2005-2010. 



VIII. APPENDIX A: WAVE FUNCTIONS AND NORMS 



In this appendix we shall derive alternative expressions of the eigenfunctions of the model 
and compute their norm. Let us start from eq. fl4.7p for the eigenfunctions of the Hamiltonian 



-{l/2-iE)q 



Co+ f dq' e-''"i'{Ba{q')-Ah{q')) 

J — oo 



Replacing a{q) and 6(g) by their Fourier transform, and using eq.( ]5.3p one finds 



dq' e-'^«'a(g'^ 



duo e*«^a(-cu) 



(8.1) 



(8.2) 



2 ■ ' 27ii uo- E 

and a similar expression for the integral of h{q). All the singular integrals appearing in this 
appendix must be understood in the Cauchy sense. Plugging the latter expressions into 
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dUD yields 



Co + ksai-E) - Abi-E)) + e"-- H ^ e- ^^^""^ " f ^'"^ 
2 y_oo ^TT? a; - 

(8.3 

Using eqs.dM]), KTj\\ and fOHD . the first term in the RHS becomes 

Co + \{Ba{-E) - Ab{~E)) = 9^±3^ = Re J^{E) (8.4) 
so that ip{x) is given by 

where A(-E) and B{E) are given by the eqs. fl4.27p and fl4.28p . The function (18.50 can also 
be expanded in the basis (13. 3p of eigenfunctions of Hq, i.e. 



\^e) = I doo Moo) 10^) (8.6) 

namely 



M^)= I duoi,E{^) (8.7) 
) V 27r 



The result is 



, . I — s , . 1 B(E)a(-uj) - A(E)b(-uj) , , 

^B(tu) = v^5(E-cj) Re J^(E) + ^ (8.8) 

V 27r? <jj- E 



which shows that the delocalized states, i.e. J^[E) ^ 0, have to be normalized in the 
distributional sense, while the localized states, i.e. J^{Em) = 0, have a norm given by 

|, . r \B{E^)a{-uj) - A{Em)b{-uj)\^ 
{'^eJ^eJ = / 7 r—^ (8.9) 

In the examples discussed throughout the paper the functions a{t), b{t) are phase factors, up 
to overall constants. Moreover, if the function a{t)b{—t) is analytic in the upper half-plane 
and vanishes when |t| oo. Re t > 0, then the S'-functions and the associated Jost function 
take a particular simple form if we allow for the existence of bound states, 

Sa,a = Sb,b = 1, Sa,b = ait) b{-t) , Sb,a = =^ ^(t) = 2 + a(t) b{-t) (8.10) 

The integration constants A, B, corresponding to a bound state, can be choosen as 

A{EJ = -B{EJ = -1 (8.11) 
37 



which differ with respect to (14.281) in an unimportant overall sign. The wave function (18. 5p 
also simplifies 

V>..W=r^:^-^-^- (8.12) 
and scalar product of two bound state wave functions becomes 

WE^,\i^E^J= 77-7- p. w , (8-13) 

J-00 - Emi){uJ - Em^) 

The analiticity of the Jost function J^{E) in the upper-half plane implies the dispersion 
relation 

r(E)^T^+ r'^^ (8.14) 

where J^oo is the value of J^{E) at E = +ioo. From this equation, and the fact that 
J^{Emi) = J-'{Em2) = 0, one can show that ipE^i i'Emi orthogonal. Furthermore, 
eq. fl8.i4p yields also a simple expression for the norm of ipEm 

H^eJ^eJ = T- 7^^^ = -Im riEm) (8.15) 

Finally, writing J^{E) as in eq. (17.380 . i.e. 

J^(E) = 2(1 + e e'™^^)) (8.16) 

where n{E) is the number of states, up to a constant, one derives that the norm of ipEm is 
proportional to the density of states at Em, 

{i^Eji^Ej=^T^n\Em) (8.17) 

A. Wave functions associated to the smooth and exact Riemann zeros 

The Mellin transforms of the boundary wave functions associated to the smooth Riemann 
zeros were given in eq. (l6.10p . Choosing = l,lp = 27r, ao = &o = we have 

a{t) = ^/2e^''^'\ b{t) = 1 (8.18) 

The wave functions (I8.12p become in this case, 

^ (8.19) 
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FIG. 10: Plot of \ipE,J for the energies Riemann zeros: 14.1347, 21.022, 25.0109 evaluated with 
eq. (l8.2ip for A = 60. The wave function are normalized using eq. (l8.17p . 

The integrals can be performed using the residue theorem obtaining 

1 , , , H(x — 1) 1 „ ,1 iEm 1 5 iEm o , X 

W,^ X = , ^ + 1 ^ ii^2 -]-, ^,-7rV 8.20) 

^J2 I ^ lEsL ^4 2 2 4 2 

where H{x — 1) = 1 if x > 1 and if < x < 1. One can show that y/xtpE^ — ^ as 
X oo, if 1 + e^*^*^-^™) = 0. In figl9]we plot the absolute values of fl8.20p for those energies 
that correspond to the three lowest Riemann zeros. Notice that the functions are very small 
in the classical forbidden region < a; < 1. The amplitude has a high frequency component 
common to the three waves plus a low frequency one that depends on the level. 
The wave functions associated to the exact Riemann zeros can be computed from eq.f l8.12p 
with a{t) and b{t) given by eq. (17.411) . We do not have an analytic expression for this 



39 



integral, however a numerical estimate can be obtained truncating (18.121) as 

Je^-k 27rz uj-E,n 
In figHO] we plot the result for the lowests Riemann zeros. The wave functions have some 
common features with those of fig. [9l but they also exhibit a random behaviour. 
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